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ABSTRACT 

We review in this work specific-heat experiments, that we have conducted on differ- 
ent hydrogen-bonded glasses during last years. Specifically, we have measured the low- 
temperature specific heat Cp for a set of glassy alcohols: normal and fully-deuterated ethanol, 
1- and 2- propanol, and glycerol. Ethanol exhibits a very interesting polymorphism present- 
ing three different solid phases at low temperature: a fully-ordered (monoclinic) crystal, an 
orientationally-disordered (cubic) crystal or 'orientational glass', and the ordinary structural 
glass. By measuring and comparing the low-temperature specific heat of the three phases, in 
the 'boson peak' range 2-10 K as well as in the tunneling-states range below IK, we are able 
to provide a quantitative confirmation that "glassy behavior" is not an exclusive property 
of amorphous solids. On the other hand, propanol is the simplest monoalcohol with two 
different stereoisomers (1- and 2-propanol), what allows us to study directly the influence of 
the spatial rearrangement of atoms on the universal properties of glasses. We have measured 
the specific heat of both isomers, finding a noteworthy quantitative difference between them. 
Finally, low-temperature specific-heat data of glassy glycerol have also been obtained. Here 
we propose a simple method based upon the soft-potential model to analyze low-temperature 
specific-heat measurements, and we use this method for a quantitative comparison of all these 
data of glassy alcohols and as a stringent test of several universal correlations and scaling 
laws suggested in the literature. In particular, we find that the interstitialcy model for the 
boson peak [A. V. Granato, Phys. Rev. Lett. 68 (1992) 974] gives a very good account of 
the temperature Tmax at which the maximum in Cp/T'^ occurs. 

PACS numbers: 65.40.+g, 63.50.+X, 61.43.Fs 
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1. Introduction 

After several decades of research on the subject, the universal properties exhibited 
by glasses at low temperatures |jl|, |^ (i.e., their vibrational excitations at low frequen- 
cies) continue to be a vivid matter of interest and debate 0. It is well known 0, ^ 
that glasses or amorphous solids have thermal properties (and also dielectric or acous- 
tic ones) very different from those of crystalline solids. Moreover, these properties 
are very similar among different families of glassy materials irrespective of either the 
type of chemical bonding or other structural details, hence the name "universal". At 
temperatures T < 1 K, the specific heat Cp depends approximately linear {Cp oc T) 
and the thermal conductivity n almost quadratically (k oc T^) on temperature, in con- 
trast to the cubic dependences observed in crystals for both properties which is well 
understood in terms of Debye's theory. Quantitatively, the specific heat of nonmetallic 
glasses is much larger and the thermal conductivity orders of magnitude smaller than 
those of dielectric crystals. At T > 1 K, Cp still deviates from the expected CDebye oc 
dependence, presenting a broad maximum in Cp/T^, in the same temperature range 
where the thermal conductivity exhibits a universal plateau. It is nowadays clear that 
this universal feature is closely related to an excess in the vibrational density of states 
g{v) over the crystalline Debye behavior, leading to a ubiquitous maximum in g{v)jv^ 
at frequencies ~ 1 THz which is known as the boson peak, a dominant feature in the 
vibrational spectra of glasses very thoroughly observed and studied by Raman and 
inelastic neutron scattering. 

About 30 years ago, Phillips and Anderson, Halperin and Varma |^ introduced 
independently the now well-known tunneling model (TM), which postulated the ubiq- 
uitous existence of atoms or small groups of atoms in amorphous solids which can 
tunnel between two configurations of very similar energy. This simple model of tunnel- 
ing states successfully explained the low-temperature properties of amorphous solids 
[0, though only for T < 1 K. However, the also rich universal behavior of glasses above 
1 K (the broad maximum in Cp/T^, the corresponding boson peak in vibrational spec- 
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tra, or the abovementioned plateau in the thermal conductivity) still were unexplained. 
Among the different approaches proposed since then to understand the general behavior 
of glasses in the whole range of low-frequency excitations, the phenomenological soft- 
potential model (SPM), which can be regarded as an extension of the TM, is one of the 
best accepted and most often considered. The SPM ||, |[ postulates the coexistence in 
glasses of acoustic phonons (crystalline-like, extended lattice vibrations) with quasilo- 
calized low-frequency (soft) modes. In the SPM, the potential of these soft modes is 
assumed to have a uniform stabilizing fourth-order term W, an energy scale which is 



the basic parameter of the model M, U, 10, 11, O, O]. In addition, each mode has 



its own first-order (asymmetry Di) and second-order (restoring force D2) terms, which 
can be either positive or negative, hence giving rise to a distribution of double-well po- 
tentials (tunneling states) and more or less harmonic single- well ones (soft vibrations). 
The parameter W marks the crossover from the tunneling-states region at the low- 
est temperatures to the soft-modes region above it. Indeed, W can be approximately 
determined either from the minimum Tmin in Cp/T^ {W ~ 1.8-2 /cB^mm) or from the 



position of the maximum Ti^^max in a k/T versus T plot: (W ~ 1 .6 k-BTi^^max) [0! HI 
Similarly to the TM, a random distribution of potentials is assumed: P{Di, D2) = Ps- 
For a more detailed description of the SPM, the reader is referred to the reviews of 



Refs. 14, 15 



In addition to these phenomenological models, other recent approaches [TO, 118 



[T^, ^ have focused on suggesting general scalings, correlations or universalities in the 
low-temperature properties of glasses which could hint at their microscopic origin. 

In order to gain understanding in this issue, we have conducted a series of measure- 
ments [pl| , |2^ , [23| , p^ , p^ , p6| of specific heat at low temperatures for a special family 
of glasses, simple alcohols such as ethanol |2^, |23|, propanol and glycerol [p5[] , 
which have low glass transition temperatures Tg (they are liquid at room temperature) 
and a molecular, hydrogen-bonded network. In this work, we propose a systematic 
method to analyze low-temperature specific-heat data, partly based upon the SPM. 
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Then, we employ this method for a critical comparison of the data obtained by us for 
the alcohols and, finally, we make use of the collected set of parameters as a test of 
several correlations and scaling laws suggested in the literature. 

2. Experimental 

During last years, we have performed specific-heat measurements on different glassy 
alcohols, by employing a calorimetric cell especially designed for samples which are 
liquid at room temperature and that allows to prepare in situ different solid phases. 
Firstly, experiments were conducted in a ^He cryostat, reaching temperatures down to 
~1.7 K. Later, we have used a very similar calorimetric cell within a ^He cryostat being 
able to measure the specific heat to about 0.5 K. 



In particular, we have studied [21, 22, 23| both normal hydrogenated and fully- 



deuterated ethanol, both of which exhibit a rich polymorphism, presenting stable 
(monoclinic) crystal, orientational glass (OG, an orientationally-disordered cubic crys- 
tal, obtained by quenching a rotationally-disordered plastic phase), and structural glass 



(amorphous) phases p7 



We have also investigated ||2^ the behavior of the next substance in the series of 
monoalcohols, propanol, which is the smallest one which has two different stereoiso- 
mers, 1- and 2-propanol, hence allowing us to study directly the effect of the spatial 
rearrangement of atoms on the low-temperature properties of glasses. 

Finally, we have also measured the specific heat and the thermal conductivity 



of glassy glycerol, probably the most widely studied glass-forming liquid. 

Further details on the experimental setup employed for the heat-capacity measure- 
ments, as well as on the experimental procedures followed to obtain and characterize 
the different solid phases are given in the corresponding references indicated above. 

3. Results 

The specific heat of several glassy alcohols (deuterated ethanol in either true glass 
or OG phases, both isomers of propanol, and glycerol) is shown in Figs. 1 and 2, all 



of them exhibiting the usual "glassy" behavior, with a quasilinear contribution at very 
low temperature (tunneling states) and a broad maximum in Cp/T^ (boson peak). 

First of all, we remark that both structural glass and OG (i.e., a crystal with 
orientational disorder) phases of ethanol show, qualitatively and even quantitatively, 
the same glassy features in the low-temperature specific heat p2|, the boson peak in 



inelastic neutron scattering [^, or at the glass transition ||22|. These results provide a 
quantitative confirmation of the fact that "glassy" behavior is not an exclusive property 
of amorphous solids, which simply lack translational disorder, but a more general 
characteristic of solids where any kind of disorder is able to soften the rigid vibrational 
spectrum of a crystalline lattice and/or to provide additional thermodynamic degrees 
of freedom which are somehow the basic ingredient of the glass state. 

It is also noteworthy that 2-propanol has a much larger specific heat above 1 K 
than 1-propanol (and than any other glassy alcohol). However, the reason for this 
difference is based on its significantly larger Debye contribution (see Fig. 1 and Table 1). 
This difference between both isomers of propanol also occurs in their crystalline states 
2^], which have been recently found to be different indeed, a monoclinic crystalline 



structure for 1-propanol and a triclinic one for 2-propanol [^. Therefore, the influence 
of the position of the hydroxyl (0H-) on the elastic constants of the hydrogen-bonded 
network seems to be very relevant. In contrast, the "excess" specific heat attributable 
to tunneling states and quasilocalized vibrations in glasses is much less affected by 
these changes in the atomic arrangement. 

In addition, we have concurrently measured the specific heat and the thermal con- 
ductivity of glassy glycerol |2^, and have used those combined data as a more reliable 
test of the SPM, which has been shown to successfully explain the specific heat and 
the thermal conductivity in a wide temperature range. In this work, we will only use 
the specific-heat data for a general comparison with other glassy alcohols. 

In Fig. 1, the specific heat of five glasses (four truly amorphous and one disordered 
crystal) is displayed below ~2.5 K in a typical Cp/T vs plot. In the simplest 
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version of the TM 0, ^ the random distribution of tunnehng states can be regarded 
as a constant density of two- level systems (TLS), which neglecting some logarithmic 
corrections gives a linear contribution to the specific heat Cp = Ctls^- Taking also 
into account the Debye contribution due to extended long-wavelength vibrations of the 
amorphous lattice, Coebye = Cr,T^, it is clear that Cp/T vs should plot linearly 
with an intercept Ctls at T = 0. Although this simple method has been traditionally 
used to determine Ctls and Cd, it poses some problems. As Fig. 2 reminds, there is an 
additional source of specific heat arising from the low-frequency vibrations responsible 
for the boson peak and the maximum in Cp/T^. This contribution is not completely 
negligible below 1 K, that brings as a consequence that many linear fits Cp = CiT + 
C3 found in the literature unavoidably provide cubic coefficients clearly exceeding 
the true Debye one, obtained from elastic measurements in those cases where they are 
available, i.e., C3 > Cd. Indeed, Ci and especially C3 depend on the chosen range 
for the linear fit. In order to solve these inconveniences, we propose a systematic 
method to analyze low-temperature specific-heat measurements. The basic point is 
to realize that the difficulties originate from the lower-energy side of the additional 
vibrations responsible for the boson peak. These have been well accounted for by the 



SPM |T^, |T5l as harmonic soft modes giving rise to a specific heat below the maximum 
Cp = CsmT^ = ^Ps^B (^^) • Therefore, a quadratic polynomial fit in Cp/T vs 
seems the most appropriate solution. The question however remains how to decide 
the temperature range to fit the data, with physical meaning. It is clear that the 
distribution of soft modes and correspondingly the specific heat can not grow Cp oc 
unlimitedly. Gil et al. |]TT]| proposed a gaussian distribution in the asymmetry of the 
soft potentials, which without any further fitting parameter allowed them to account 
for the specific heat, thermal conductivity and vibrational density of states g{v)/v'^ in 
the whole relevant range, including the "boson peak". At least for practical reasons, 
let us assume that distribution function. It is easy to find that the simple CsmT^ 
approximation starts deviating ~ 5% from the Cp{T) curve accounting for experimental 
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data approximately at T > 0.75W ~ ^T^in- Therefore, we suggest to fit specific- 
heat data in a Cp/T vs representation by using a quadratic polynomial Cp = 
CtlsT + CdT^ + CsmT^ in the temperature range < T < |Tmin. The results 
from these fits for the studied glassy alcohols are shown in Table 1. As a proof of 
consistency, we want to mention that the obtained Cd for glycerol agrees better than 



1% (hence within the experimental error) with the Debye coefficient estimated p5| 
from elastic measurements. Unfortunately, this is the only glass from those studied 
here for which elastic data are available. Nevertheless, we believe that the method 
proposed is a reasonable alternative to determine the Debye coefficient of a glass from 
low-temperature specific-heat measurements, especially when elastic data are lacking. 
The so-obtained values of Cd and 9d have been used in Fig. 2 to scale Cp data for the 
reasons discussed below. 

For sake of completeness, we also show in Table 1 the SPM parameters W (deter- 
mined from W ~ 1.8 k-QTmin [0, 15 1, with the exception of glycerol which has a flat 



minimum and has been better determined from available thermal-conductivity data 



2^]) and Pg (determined from the given expression for Csm)- If is to be stressed that 
all studied glasses present comparable values of Ps (basically, the distribution density 
of quasilocalized excitations, either tunneling states or soft modes), not existing any 
significant difference between amorphous glasses and orientational glasses (OG). 

4. Discussion 

In this section, we will make use of the compiled set of data for the studied glassy 
alcohols taken as a model system to critically discuss several correlations or scaling laws 
which have been suggested in the literature to be universal for glasses or amorphous 
solids. 

First, we will address very briefly the scaling proposed by Liu and Lohneysen [|T^ . 
They suggested a general correlation between the mechanisms leading to the Cp/T^ 
maximum in crystalline and amorphous solids. They plotted the height Pc of the max- 
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imum in Cp/T^ (i.e., {Cp/T^)raax in our notation) versus its position T^ax for a wide 
set of materials, mainly amorphous polymers and metals, as well as typical network 
glasses, together with their corresponding crystalline solids. They found an approxi- 
mate general correlation oc T^l'^, somehow indicating a close relation between the 
Cp/T^ maxima in glasses and crystals. However, we have included in such a graph 
(see Fig. 2 of Ref. ||2^) specific-heat data for molecular glasses and crystals (both 
our hydrogen-bonded materials and van der Waals glasses from Lindqvist et al. ||30||), 
finding that these molecular solids deviate systematically about a factor 5 from the 
general scaling proposed by those authors. Moreover, for the different phases of H- and 
D-ethanol we have plotted {Cp/T^) / Pc vs T/T^ax, as also suggested by Liu and von 



Lohneysen |T6[. We found (see Fig. 11 of Ref. [Q) that data for the disordered solids 
(either glass or OG) superimpose in the whole temperature region but, in contrast, the 
curves for the ordered crystals showed a far narrower shape. This fact points out again 
to the different nature of the low-energy vibrational spectra of glasses and crystalline 
solids. 

Other authors have suggested that the low-energy excitations of glasses may be 
correlated with the fragility of the glass-forming liquid, a parameter proposed by An- 



gell |32| to characterize how fragile or strong is a supercooled liquid to resist the 
structural degradation induced by temperature, and which is usually measured by the 



deviation of the shear viscosity from an Arrhenius law. Zhu [jT9[ has suggested a general 
correlation between the density of tunneling states and the fragility, presenting data 
for a variety of glasses. More fragile glasses would have a larger number of minima on 
the potential-energy hypersurface, hence explaining a higher density of tunneling states 



and a higher value of Ctls- On fhe other hand, Sokolov et al. have correlated 
the fragility of the system with the excess of vibrational excitations normalized to the 
Debye level. Specifically, they showed that the ratio (Cexc/C'D)max = (C'p/Co — l)max 
decreases with increasing fragility. More recently, Zhu and Chen |]19| have argued 
against the universality of such a correlation. 
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Unfortunately, accurate fragility values using the same criteria are not available 
for the four glass-forming liquids studied in this work. Nevertheless, it can be seen 
from viscosity data |]3ll , ^ that all alcohols possess an intermediate degree of fragility. 
Liquid ethanol seems to be slightly more fragile than glycerol, and 1-propanol is less 
fragile. There are no published data for 2-propanol to our knowledge. In this frame- 
work, the large difference in Ctls between ethanol and glycerol (almost one order of 
magnitude, see Table 1) for two systems of very similar fragility claims against the 



correlation proposed by Zhu [0. Moreover, 1-propanol has a contribution to Cp from 
tunneling states almost a factor of 3 higher than that of glycerol, being less fragile. The 
correlation with (Cexc/C*D)max is more difficult to ascertain. The obtained values are 
all around {Cexc/Cj))^.^^ ~ 0.5 as expected for intermediate fragilities, but the smaller 
deviations do not show up the expected trend. Neither suggest normalized Cp/Co 
curves shown in Fig. 2 any kind of universal behavior. Furthermore, the very similar 
values of Ctls and (Cexc/CD)max for either true glasses (quenching a supercooled liquid) 
or orientational glasses (quenching a plastic crystal) of ethanol cast doubts about the 
relevance of the fragility of the glass-forming liquid to understand the low-temperature 
properties of glasses. 

Finally, we would like to consider the interstitialcy model proposed by Granato 



33, 34 1 . According to this model, liquids can be considered as crystals containing 



a few interstitials in thermal equilibrium, which become frozen in the glassy state. 
Self-interstitial resonance modes would be the physical realization of the soft modes 
and tunneling states of the SPM, giving rise to the boson peak. For the simplest 
approximation, taken a single frequency for the resonance modes, the maximum in 
CpjT^ is predicted to appear at a temperature T^ax ~ 6*0/35. Several authors have 
tested this relation for different collections of experimental data, finding empirically 
^max ~ 6*0/40 ||3^ or Tmax ~ 6*0/38 It cau be seen in Table 1 that this correlation 
works very well for all glassy alcohols studied by us. On the contrary, the maxima of the 



corresponding crystalline states p2| , p4| , p5[] take place at higher reduced temperatures 
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(^max/^D ~ 1/25 for these alcohols), as it seems always be the case (see Fig. 3.3 
in Ref. M). However, we have tested this correlation between T^^x and 6*0/35 from 



published data |^ for another model system: vitreous B2O3 submitted to different 



thermal treatments producing significant changes in mass density, elastic constants. 



etc. Variations of up to 17% in Tmax and 6b were achieved with this method. We 
find that, although (Cexc/CD)max ranges uncorrelatedly from 1.18 to 1.70 for the seven 
different B2O3 glasses, the ratio between T^ax and remains constant for all of them: 
O-o/Trnax = 51 ± 1. Let us notice that for other oxide glass as Si02, 6-D/Tma.x = 46 
[pO|] . So, it may well be that the exact T^^x/Od ratio could somewhat depend on the 
kind of glass, what is not in conflict with the interstitialcy model |Q, but for a given 
"system" , fixing some secondary parameters, the boson peak and Debye temperatures 
seems to be clearly correlated. 

5. Conclusions 

In summary, we have reviewed and comparatively discussed specific-heat experi- 
ments, that we have conducted on different hydrogen-bonded glasses during last years: 
normal and fuUy-deuterated ethanol, 1- and 2- propanol, and glycerol. We have pro- 
posed a systematic method partly based upon the soft-potential model to analyze these 
low-temperature specific-heat measurements, and to test several universal correlations 
and scaling laws suggested in the literature. In particular, we have found that the 
correlation between the temperature r^ax at which the maximum in Cp/T^ occurs and 
the Debye temperature proposed by the interstitialcy model for the boson peak is 
very well fulfilled. 
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FIGURE CAPTIONS 



FIGURE 1: Low-temperature specific heat plotted as Cp/T vs for several glassy 
alcohols. Solid hnes are fits to quadratic polynomials (see text). 

FIGURE 2: Low-temperature specific heat of the same glasses in Fig. 1, scaled to 
the cubic Debye contribution CDebye; ^ ^ function of temperature normalized to O-q. 
The arrow indicates the position of the maximum at T = ^d/35 predicted within the 
interstitialcy model. 
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TABLES 





H-ethanol 


D-ethanol 


-1 1 

1-propanol 


r\ 1 

2-propanol 


1 1 

glycerol 




glass 




1 

glass 


OG 


glass 


glass 


glass 


Pmoi (g/mol) 




46.1 


52.1 




60.1 


60.1 


92.1 


(K) 


95 


95 


95 


95 


98 


115 


185 


T^n (K) 


2.3 


2.6 


2.1 


2.3 


1.8 


1.6 


2.0 


Tmax (K) 


6.1 


6.8 


6.1 


6.4 


6.7 


5.0 


8.7 


\S^pl ^ J max 


2.4 


2.2 


2.8 


2.6 


2.7 


3.6 


1.4 


(mJ/molK'^) 
















Ctls (mJ/molK^) 
Cd (mJ/molK^) 


1.2 


1.27 


1.05 


1.13 


0.424 


0.516 


0.157 


1.55 


1.45 


1.80 


1.72 


1.77 


2.54 


0.855 


(mJ/molK6) 


0.0432 


0.0288 


0.0572 


0.0419 


0.0367 


0.0845 


0.0139 


W (K) 


4.1 


4.7 


3.8 


4.1 


3.3 


2.9 


4.3"f 


(mol-^) 


4.0x10^9 5.2x10^9 


3.6x10^9 3 


.8x10^9 


1.1x10^9 


1.4x10^9 


1.6x10^9 


{C^xc/ C'D)max 


0.55 


0.52 


0.56 


0.51 


0.53 


0.42 


0.64 


Ojy (K) 


224 


229 


213 


217 


236 


209 


317 




37 


34 


35 


34 


35 


42 


36 



Table 1: 

Measured data and fit parameters obtained for several studied glassy alcohols (see 
text). * Data for the glass phase of H-ethanol were taken in a ^He-cryostat, only 
down to 1.7 K. ^ The value of W for glycerol has been determined from thermal 



conductivity data [g5 
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Ramos et al., FIG 2. 




